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弹性 支撑 功能 梯度 压 电 多 了 筷 微 加 柱 元 的 
自由 振动 分 析 


刘 文 光 ! ,庞大 !, 吕 志 鹏 !, 刘 超 *, 张 宇航 


(1. 南昌 航空 大 学 航空 制造 工程 学 院 ,330063 南昌 ;2. 哈尔滨 工业 大 学 (深圳 ) 理 学 院 ,150001 深圳 ) 


摘 要 : 旨 在 研究 热 - 力 - 电 载荷 下 弹性 支撑 功能 梯度 压 电 多 有 孔 微 圆 柱 党 的 自由 振动 。 首 先 ,建立 弹 
性 支撑 功能 梯度 压 电 多 孔 微 圆柱 党 动力 学 模型 ;然后 ,应 用 三 阶 剪 切 变形 党 体 理论 和 修正 的 偶 应 力 
理论 ,推出 弹性 支撑 功能 梯度 压 电 多 有 孔 微 圆柱 党 模 态 频率 的 解析 解 ; 最 后 ,通过 数值 算 例 分 析 了 微 
柱 党 模 态 频率 的 影响 因素 。 结 果 表 明 :Pasternak 弹性 支撑 比 Winkler 弹性 支撑 更 有 利于 提高 微 圆 
柱 党 的 模 态 频率 ;改变 弹性 支撑 的 刚度 系数 、 轴 向 力 、 外 加 电压 、 孔 隙 分 布 、 材 料 体积 分 数 指 数 和 结 
构 尺 寸 可 调节 微 圆 柱 党 的 模 态 频率 ;孔隙 体积 分 数 越 大 ,温度 或 轴 向 力 对 模 态 频率 的 影响 越 大 ,而 
电压 对 模 态 频率 的 影响 则 越 小 ;不 同 材 料 指 数 下 , 增 大 孔隙 体积 分 数 对 模 态 频率 的 影响 趋势 不 同 ; 
弹性 支撑 会 减弱 温度 、 轴 向 力 和 电压 对 模 态 频率 的 影响 ,对 薄 圆 柱 党 或 短 圆 柱 党 模 态 频率 的 影响 较 
为 显著 。 
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Free vibration analysis of a functionally graded piezoelectric 
porous cylindrical microshell resting on elastic foundation 


LIU Wenguang! , PANG Lei ,LU Zhipeng ,LIU Chao? ,ZHANG Yuhang! 


(1. School of Aeronautical Manufacturing Engineering , Nanchang Hangkong University ,330063 Nanchang , China ; 


2. School of Science , Harbin Institute of Technology , 150001 Shenzhen, China) 


Abstract : This article studies the free vibration of a functionally graded piezoelectric porous cylindrical 
microshell resting on elastic foundation under thermo-electro-mechanical loads. First, the dynamic model of 
the functionally graded piezoelectric porous cylindrical microshell resting on elastic foundation was 
established. Then the analytical solution of modal frequencies of the functionally graded piezoelectric 
porous cylindrical microshell resting on elastic foundation was derived using the third-order shear 
deformation shell theory and the modified couple stress theory. Last the factors affecting model frequency 
were discussed by numerical case studies. The results show that the Pasternak foundation is preferred to 


improve the modal frequencies of the shell to the Winkler foundation. The modal frequencies of the 


收 稿 日 期 : 
基金 项 目 : 
通信 作者 : 
引用 格式 : 


2021-07-01 修 回 日 期 :2022-11-04 

国家 自然 科学 基金 资助 项 目 ( No. 51965042) 

刘 文 光 ,教授 。E-mail:liuwgl4@ nchu. edu. cn 

刘 文 光 , 庞 舌 , 吕 志 鹏 ,等 弹性 支撑 功能 梯度 压 电 多 孔 微 圆柱 壳 的 自由 振动 分 析 [ 相 . 应 用 力学 学 报 ,2024,41(2) :411-421. 

LIU Wenguang,PANG Lei,LU Zhipeng, et al. Free vibration analysis of a functionally graded piezoelectric porous cylindrical microshell resting 
on elastic foundation| J ] . Chinese journal of applied mechanics ,2024 ,41(2) :411-421. 


412 应 


力学 和 学报 


P41 È 


cylindrical microshell could be adjusted via changing the stiffness coefficients of the elastic foundation, 


axial load, applied voltage, porosity distribution, material volume fraction indexes, and structure 


dimensions. The larger the porosity volume fraction, the greater the influence of temperature and axial load 


on the modal frequency, and the smaller the influence of applied voltage on the modal frequency. The 


change trends of modal frequency with increasing porosity volume fraction are obviously different with 


different material volume fraction indexes. The elastic foundation can weaken the influence of temperature, 


axial load and applied voltage on the modal frequency , and its influence on the modal frequency of thin or 


short cylindrical shell is more significant. 


Key words functionally graded piezoelectric material; elastic foundation; porosity ; modified couple stress 


theory ; higher-order shear deformation theory 


压 电 材料 因为 呈 脆 性 ,一 般 与 弹性 基体 结合 组 
成 多 层 结构 。 然 而 ,多 层 结构 的 分 界面 会 导致 应 力 
集中 , 易 产 生 蠕 变 、 裂 纹 和 脱 层 等 各 种 失效 六 。 功 能 
梯度 材料 (functionally graded materials, FGMs ) K| -E 
有 和 承受 高 温 、 去 除 界面 失 配 等 优点 ,能 克服 以 上 一 系 
IRS TEWAK FGMs 引入 到 压 电 器 件 的 
设计 中 ,提出 了 功能 梯度 压 电 材料 (functionally 
graded piezoelectric materials, FGPMs) 。 因 为 FGPMs 
具有 高 机 电 系 数 、 高 柔性 度 和 高 届 服 强度 的 优点 ,在 
压 电 振 动能 量 俘获 等 方面 具有 很 多 优势 。 
围绕 FGPMs 结构 的 振动 、 主 动 控制 和 届 曲 等 问 
题 ,研究 者 做 了 很 多 工作 。 如 果 结 构 尺 寸 小 至 微 纳 
米 级 时 ,其 力学 性 能 与 宏观 结构 相 比 具有 显著 差 
异 中 。 研 究 者 针对 这 一 问题 提出 了 描述 微观 结构 力 
学 性 能 的 修正 偶 应 力 理论 '" 。 基 于 修正 的 偶 应 力 理 
论 ,SIMSEK 等 ”探究 了 含 孔隙 FGMs 微 板 在 移动 载 
荷 作用 下 的 受 迫 振动 ,KHORASANI 等 ”运用 无 网 
格局 部 Petrov-Galerkin 法 分 析 了 FGMs f B) 8$ 25 
曲 行为 ,AMIR 等 H6 T Pasternak 弹性 支撑 微米 
梁 的 热 致 自由 振动 。 考 虑 大 搁 度 行为 以 及 小 尺度 效 
应 ,KOMIJANI 等 "建立 了 FCPMs 梁 在 多 场 耦合 载 
答 作 用 下 的 非 线 性 时 变动 力学 模型 。 结 合 高 阶 剪 切 
变形 壳 体 和 修正 的 偶 应 力 理论 ,DEHSARAJI A 
研究 了 考虑 厚度 伸展 效应 的 FGPMs 微 /纳米 圆柱 过 
在 均匀 温度 场 和 电场 作用 下 的 自由 振动 响应 。 

因 材 料 制备 或 者 加 工 制造 等 原因 ,FGPMs 结构 
内 部 不 可 避免 地 会 出 现 纳米 孔洞 或 孔隙 。 为 了 
更 加 准确 地 分 析 结 构 的 振动 行为 ,有 必要 在 动力 学 
建 模 中 考虑 孔隙 对 FGPMs 结构 振动 的 影响 。 考 虑 
三 项 惯性 力 并 应 用 Kirchhoff 板 理论 ,WANG ”分 析 
了 平移 状态 下 含 均匀 分 布 孔 际 FGPMs 板 的 非 线 性 
振动 啊 应 。 采 用 改进 的 四 变量 前 切 变形 板 理 论 ， 


BARATI 等 .研究 了 多 种 边界 条 件 下 FGPMs 多 孔 
板 的 热 - 力 - 电 振 动 特性 。 应 用 非 局 部 弹性 理论 和 高 
阶 剪 切 变 形 板 理论 ,ZENKOUR 等 “探究 了 热电 载 
Ta ERI TEILE FGPMs 纳米 板 的 屈曲 行为 。 考 虑 
均匀 和 非 均 匀 孔 阶 分 布 ,SU 4&9 用 改进 的 变 分 法 
分 析 了 弹性 支撑 FGPMs 梁 的 自由 振动 和 瞬 态 响应 。 
此 外 ,考虑 结构 与 周围 弹性 介质 间 的 作用 ,研究 者 先 
后 提出 了 模拟 弹性 介质 的 Winkler 弹性 支撑 和 
Pasternak 弹性 支撑 模型 '" 。 基 于 弹性 支撑 模型 , 陆 
续 有 研究 者 讨论 了 弹性 支撑 对 结构 振动 特性 的 
影响 1。 

尽管 研究 者 对 含 孔 隙 FGPMs 结构 的 振动 问题 
开展 了 大 量 的 研究 ,但 是 鲜 有 文献 以 修正 的 偶 应 力 
理论 为 基础 ,探讨 弹性 支撑 FGPMs 多 孔 微 圆 柱 壳 的 
振动 啊 应 。 本 研究 以 Pasternak 弹性 支撑 FGPMs 多 
孔 微 圆 柱 壳 为 对 象 ,基于 三 阶 剪 切 变形 理论 建立 模 
型 的 几何 方程 以 修正 的 偶 应 力 理论 反映 微 纳米 结构 
所 产生 的 尺度 效应 ,进而 推导 出 FGPMs 微 圆 柱 壳 的 
多 场 耦 合 致 振动 方程 ,对 FGPMs 微 圆 柱 壳 的 自由 振 
动 特性 进行 分 析 。 研 究 结果 可 为 压 电 微 纳米 器 件 的 
动力 学 设计 提供 一 定 的 支撑 。 


S 


1 FGPMs 微 圆柱 壳 模 型 


如 图 1 所 示 的 Pasternak 弹性 支撑 FGPMs Z 4L 
微 圆柱 壳 , 壳 长 为 L, 壁 厚 为 h, 中 面 半 径 为 R。 弹 性 
支撑 均匀 分 布 在 微 圆 柱 壳 的 外 表面 ,其 中 包括 弹性 
刚度 系数 为 瓜 的 弹簧 层 和 剪 切 刚 度 系数 为 户 的 剪 
切 层 。 当 剪 切 层 刚 度 系数 所 等 于 0 H}, Pasternak 弹 
性 支撑 转化 为 Winkler 弹性 支撑 。 在 圆柱 壳 的 中 面 
上 建立 坐标 系 ,x、9、z RI y 轴 分 别 表 示 壳 体 的 轴 向 、 
环 向 、 径 向 和 中 性 层 切线 方向 ,ww 和 分 别 表 示 中 
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面 上 任意 点 沿 x\y 和 = 方向 的 位 移 。 


FGPMs- II 


FGPMs- 
图 1 弹性 支撑 功能 梯度 压 电 多 孔 微 圆柱 元 模 型 


Fig.1 Model of an FGPM porous cylindrical microshell 


resting on elastic foundation 
fr [PEE E AERE Hs E BISHER JE 28 Us PALA TR 
为 PZT-4, 外 表面 为 PZT-5H, 材 料 特性 如 表 1 所 示 。 
假设 材料 属性 沿 厚度 方向 服从 竹 律 分 布 ,同时 考虑 
均匀 分 布 孔 际 和 非 均 匀 分 布 孔 际 两 种 孔隙 分布, 则 
微 圆柱 壳 的 有 效 材料 属性 表达 式 为 


Pi (P, -P)I" +P- AP, P) (1) 


PyG) = C, - POT" +P, - P, Po (1 - 2181) 


(2) 
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式 中 :P; 和 Pi 分 别 为 FGPMs- 了 型 均匀 分 布 孔 隙 和 
FGPMs- 工 型 非 均匀 分 布 孔隙 下 微 圆柱 壳 的 有 效 材 
料 属性 ;P; 和 P, 分 别 为 壳 体 内 外 表面 的 材料 属性 ; 
T =z/h +1/2;N 为 PZT-4 的 体积 分 数 指数 ;ev 为 孔 
隙 体积 分 数 。 
假设 微 圆柱 壳 受 到 热 - 力 - 电 载 荷 的 作用 。 如 网 
1 所 示 ,在 微 圆柱 过 内 外 表面 施加 热 载荷 Ti; 沿 轴 向 
施加 静态 轴 向 力 IN, IN, D AE S ALI S, IN, DA f 
表示 轴 向 受 压 ;在 内 外 表面 间 施 加 电 载 荷 Bu ,Gu 为 
正 表 示 电 势 方向 与 微 圆 柱 壳 极 化 方向 相同 ,Gu 为 负 
则 表示 相反 。 
假定 温度 变化 仅 出 现在 壳 体 厚度 方向 ,根据 一 
维 傅 里 叶 热传导 方程 ,在 壳 体 厚度 方向 上 的 温度 梯 
ET URRH 
+ + | (3) 


T(z) = T, -7,( 
式 中 ,7,,=7 -TTA T, 4) AENA E o 
假设 圆柱 壳 材 料 沿 厚度 方向 由 内 向 外 表面 极 
化 ,电热 沿 厚度 方向 呈 二 次 分 布 ,满足 Maxwell 电学 
J m 


zo 
Q,y,z,D- ! 


*té(z)e,y,D (4) 


AP: Ter 性 变形 而 产生 的 感 生 
BR; é (2) = - (2) 。 


表 1 PZT-4 和 PZT-5H kija pe o 
Tab.l Material properties of the PZT-4 and PZT-5H ''??* 
属性 PZT-4 PZT-5H 
弹性 常数 /GPa & =139 c 577. 8,64 =74, c3 0115,64 525.6, ey 5126,05 579. 1,65 583.9, 


cq, 230.6 
压 电 常 数 /(C. m7?) 
frt it EU (107? F » m7!) 


SUSERE/ (100 Nm ?-K^!) 


热 释 电 /(10-C.m-?.K-!) Pp1=5.4,p3 =5 

密度 /(kg * m7?) p =7 600 

2 FGPMs 微 圆柱 过 运动 方程 
基于 三 阶 剪 切 TRIE ,圆柱 壳 的 位 移 

场 可 表示 为 
Az 
u, =u +z,- i 
u, =v +z, - PPS 
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ej; =12.7,e3 = -5.2,e3 =15.1 
sj 26.46,54, =5.62 


Bii =5.19,B3 24.8 


ca 2117 ,c44 223 ‚ce 223.5 
eis =17,e3 = —6.5,e4, 223.3 
sı =15.05, s33 213.02 
Bi, =5. 19 „B33 24.8 
pı =5.4,p3 =5 


p =7 500 


(5) 
式 中 :下 标 中 的 逗号 表示 偏 导 ;ui 、ws 和 ww 分别 为 党 
内 任意 点 沿 x、y 和 zz 方向 的 位 移 ;yw, 和 ,分 别 表 示 
中 面 法 线 沿 x 和 7 方向 的 转角 六 表示 时 间 。 
假设 圆柱 党 处 于 平面 应 力 状态 ,本 构 关系 为 
o = Ce - eE - ATB (6) 
D = e'g +sE + ATp (7) 
式 中 :A7=7(z) -7,, 假 定 圆 柱 过 的 外 表面 温度 7, 
为 零 热 应 力 状态 (7, =300 K);o e E,D 和 分 别 
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表示 应 力 ,应变 .电场 .电位 移 和 热 释 电 向 量 ;C、.eB 
FU s 分 别 表示 材料 弹性 常数 . 压 电 常数 、 热 模 量 和 介 


电 常 数 矩 阵 。 
由 修正 的 偶 应 力 理论 得 圆柱 壳 应 变 能 U 为 
L 27R .h/2 


jade je M: $ 
72 JoJo Jan O9 + Oey + TuYs 


Ti Yi + Ty Yy + MXi tm + MX + 
"noc + MX bay DLE - 
D,E,- D.E.)dxdydz (8) 
式 中 :cc 和 分别 为 微 圆柱 壳 x 和 > 方向 的 正 应 
JITA xy Wi VES BRI 7J sr, IE m, 29 z 2r ILES] D E 
Jd ie, 6, FI y, 、y,、y: 分 别 为 相应 方向 上 的 正 应 变 
BW ZEE, ERI E.g ER A SD, D, A D, HE 
DIEA HE ue mij =*,y,z) 分 别 为 对 称 旋转 梯 
度 分 量 和 高 阶 应 力 分 量 , 且 有 
m; = 2l eus (9) 
式 中 :1 为 描述 小 尺度 效应 的 材料 特征 长 度 ;co6. 为 微 
圆柱 壳 的 有 效 弹 性 常数 , 且 有 Cese = cu。 
微 圆 柱 壳 的 动能 K 为 


L 2mR ,.h/2 
| pC) Gi id + ipdedydz (10) 
-h/ 


Beo. 


式 中 ,p(z) 为 微 圆柱 壳 的 有 效 密度 。 
热 载荷 和 电 载荷 做 的 功 本 为 
w= -Ef f Ns ram 
(Ns) (+ Ng) w^, ] dxdy (11) 
式 中 :Ne 和 Na 分别 为 温度 在 > 和 方向 引起 的 热 
载荷 ;Ne 和 Ne 分 别 为 外 加 电压 BTE x 和 7 方向 引 
起 的 电 载荷 。 
静态 轴 向 力 做 的 功 WW 为 
W, =- | o +w? ddrdy — (12) 
忽略 弹性 支撑 的 阻尼 和 惯性 力 ,弹性 载荷 做 的 
D WH 


ES 2 2 2 
Wss- 3 f f, Teu? +k Cu, e wldrdy (13) 


因此 , 微 圆 柱 壳 的 能 量 函 数 可 表示 为 
II K-U-W,-W, +W; (14) 
通过 Hamilton 公式 可 推出 圆柱 壳 运 动 方程 


E o , 
Lo u > * (1 E si. E Na. T N 十 


x "ek h? 
i En ) (15) 
2 Ys T Y, F 2R P zu zy 
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Af, AL- 
hv 3h y (1 -jy T SAI + Noya v 
1 Q, 1 
pe Ys 十 了 * R + igl- T N,v 。 x2 
1 
PPM Ya. ul Tu ij Ya) ( 16) 
AL. ü 4I 164. . » " 
ap s F vy) 十 人 = a] Paa " p, I, w 一 
jh do oer. 
gpi Chas Pup ag Bars bab. poA 


Ny 4 
RT D ue e roam Ep + Nw 十 


1 
A Y.) Y -Y +Y, 一 了 2) + 


2R Xz,y ynx 2 YY, XY XX, Xy XY ,XX Xy.yy 


4 
pa Hes * H2 t k,Cw ,, T wy? -kuw + 


ET 

h? P eis S P ous t T ss Ec s 

Noi t Np) W nx x CN, T Np) W yy (17) 
8 16, \; 16 4 E 

(i = zp“ rd + Pa - pl), + 
AT, Y. 4 

U kn SIL = M... 十 M, ， x nu a 


4 1 
Q. + n^ bod etd re LILY,. + 


ELA zy 


2 
p2 Pes P a des angl Puch Sayy) t 
z ha ) (18) 
2 Hg F Hay E ah? Lx * Toss 

8 16,; AL, Y. 
(5 Tay yh] + (1 up t 

16 4 ，\. 2 
[ox "s ah]. = an T hd 
4 4 

Ms t M a z ag dau 十 ou) E Q, T n t 
1 1 4 
3€ Y. x Yy S Jos dl FR PE nu d 
La yag ) (19) 
2 HA, + H. 。 F h? BS u Bia s Pra 


h/2 


f TEG) D., e £()D,, -22D.1dz =0 (20) 


各 经 典 和 非 经 典 的 力 和 力矩 ,以 及 广义 惯性 党 
数 的 具体 表达 式 为 
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h/2 


N,N. Np =f Cg, 0, ,0,,) dz 
xx 25 yy 3^ xy ET, xx 3€ yy 5 U xy 


h/2 


M., M, M, = | 1L Urso, o, Didi 


h/2 


CSSS) =f Cr ar Ode (21) 
Ex 3 7 yy 9 7 ay M2 xx 29V yy? V xy 


SENE en „p Tas T a) dz 
(了 ,了 , Y, Y, ,7 Y 


ax 7 yyt azz 7 oxy?7 oxi 360 yr 


h/2 


um | "ML I, s Maz IB. sM yz m, dz 
-h 


h/2 


CHa, HO i MAC 


xz? 


m, )zdz 
| (22) 
PP PP PP) 


Xx Uo yy? omi  xy?!7 xz 7 yz 


h/2 


= Í ( m m m m m m )z dz 
xx 9 Uy 9 zz p HU y s PIU gg 9 "yz 
-h/2 » ) ? 
h/2 


Gm yer Gn, m, z dz 
-h/2 


xz 35 yz 


h/2 


Id ub bull = Í ; OO uz g pg dz 
一 AL 


(23) 
3 FGPMs 微 圆柱 壳 模 态 频率 求解 
假设 微 圆柱 壳 的 两 端 简 支 ,那么 在 任意 单一 模 


态 下 运动 控制 方程 (15) ~ (20) 的 解 可 假设 为 


u(x,y,t) = ue cos( A,,x2c08Cny/R) , 


UC x,y,D = v," sinCA, x)sin( ny/R) , 

WCx,y,D = w,, e" sinCA, x) cos ny/R), 

V.Ox,y,D = Wanne cos CA, x) cosl ny/R) , 

Wx,y,t) = Ure" sint A, sinCny/R) , 

QO,y,D- Q9, € " sin CA x) cosC ny/R) (24) 
式 中 :A mm L;m 和 分别 为 微 圆柱 党 的 轴 疝 半 
波 数 和 环 向 波 数 ;ww sous >W mn as Wynn 和 gi 为 振 
型 系数 ;w 为 角 频 率 。 

将 方程 (24) 代 入 运动 控制 方程 ,可 推出 微 圆柱 
壳 的 模 态 频率 方程 

(K - T,K, -NK, -DK -ow M)d =0(25) 
式 中 :K 为 刚度 矩阵 ;Ki、K， 和 KK, 分 别 为 温度 梯度 
7 静态 轴 向 力 N, 和 外 加 电压 B。 对 应 的 系数 矩阵 ; 
M 为 质量 和 矩阵 ;4 为 振 型 系数 向 量 。 

为 便于 分 析 讨 论 ,将 角 频 率 .静态 轴 向 力 和 弹性 
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支撑 的 刚度 系数 量 纲 归 一 化 , 即 
N = wR /p/E 
_R V3 7 vN, 


N, : 
EN? 
— DER (1 - v)? (26) 
ui Eh? 
DER (1-v) 
dii Eh? 


xXh E. p 和 分 别 为 PZT-4 ARPER EE EANA 
松 比 。 

为 验证 模型 的 准确 性 ,首先 将 模型 退化 为 纯 
BaTi0; 压 电 圆柱 党 。 图 2 比较 了 应 用 本 研究 方法 和 
应 用 文献 方法 所 得 的 结果 。 与 文献 不 同 的 是 ,本 研 
究 模 型 应 用 了 三 阶 剪 切 变形 理论 ,并 通过 修正 的 侦 
应 力 理论 反映 了 微 纳米 结构 所 产生 的 尺度 效应 。 分 
析 结 果 发 现 ,本 研究 模型 所 得 结果 与 文献 [26 ] 和 
[27] 结 果 误 差 在 2% 左右 , 且 最 大 误差 小 于 4% , 验 
证 了 本 研究 所 采用 方法 的 合理 性 。 进 一 步 将 现 有 模 
型 简化 为 纳米 圆柱 壳 , 表 2 结果 与 文献 对 比 表 明 本 
文 模型 所 得 结果 与 文献 [11 ] 和 | 24 | 结果 比较 吻合 。 

2 000 


1 800 
1600 


一 一 文献 26], R/h=50 
一 人 一 文献 [26], R/h-100 
一 e 一 文献 [26], R/h=200 
一 一 文献 [27], R/h=50 


、1400| ”一 一 文献 [27], R/h=100 
加 1 200| 一 一 文献 [27], R/h=200 
m 一 上 本 研究 , RIh=50 
gx 1000| 一 + 一 本 研究 , R/h=100 
I$ soo 一 一 本 研究 , R/h=200 


L/R=2, m-1 


0 2 4 6 8 10 12 14 16 18 20 
环 向 波 数 n 


图 2 压 电 圆柱 党 模 态 频率 的 比较 


Fig.2 Comparison of modal frequencies of 


piezoelectric cylindrical shells 
4 “分 析 与 讨论 


如 无 特殊 说 明 , 本 研究 微 圆 柱 壳 的 几何 尺寸 分 
别 取 :R=1 um,L/R =1,h/R =0.02; 材 料 特征 长 度 
1=3h; 了 和 孔隙 体积 分 数 e。=0.5 ,体积 分 数 指数 N=1; 
外 加 电压 Bu = 1 V ,无 量 纲 静 态 轴 向 力 N=0.5, 热 
载荷 T, = 300 K; 无 量 纲 刚 度 系 数 (K,,K,) = 
(1000,1), 
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表 2 不 同 厚度 下 纳米 圆柱 壳 无 量 纲 模 态 频率 的 比较 
Tab.2 Comparison of dimensionless modal frequencies of 


cylindrical nanoshells under different thicknesses 


l=0 l=h 
h/R n 


H 力学 学 报 


文献 [11] 文献 [28] 本 研究 ”文献 [11] 文献 [28] 本 研究 


0.1 1 0.978 0.932 0.936 1.179 1.135 1.163 
2 0.820 0.775 0.781 1.144 1.085 1.164 
3 0.746 0.711 0.717 1.308 1.241 1.364 
0.2 1 1.085 1.043 1.051 1.645 1.608 1.689 
2 1.009 0.967 0.974 1.771 1.690 1.908 
3 1.090 1.051 1.055 2.201 2.086 2.439 
0.3 1 1.211 1.173 1.183 2.057 2.045 2.210 
2 1.206 1.160 1.167 2.242 2.159 2.566 
3 1.388 1.337 1.334 2.785 2.652 3.271 


4.1 环 向 波 数 的 影响 
表 3 研究 了 环 向 波 数 对 微 圆柱 壳 无 量 纲 模 态 频 
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率 的 影响 。 可 以 得 出 :在 低 环 向 波 数 下 ,FGPMs- T 
型 微 圆柱 壳 的 模 态 频率 大 于 FGPMs- 荆 型 微 圆柱 壳 
的 模 态 频率 ,而 环 向 波 数 较 大 时 则 相反 。 表 明和 孔 际 
的 分 布 形 式 对 环 向 振动 的 影响 较 大 。 无 论 是 
FGPMs- 了 型 还 是 FGPMs-I 型 ,无 弹性 支撑 和 
Winkler 弹性 支撑 下 圆柱 壳 的 基 频 均 在 (m,n) = (1, 
3) 处 ;而 Pasternak 弹性 支撑 下 圆柱 党 的 基 频 出 现在 
(m,n) =(1,1) 处 。 究 其 原因 , Winkler 弹性 支撑 只 
影响 微 圆柱 完 的 径 癌 振动 ,而 Pasternak 弹性 支撑 不 
但 影响 微 圆柱 壳 的 径 癌 振动 ,还 影响 环 向 振动 。 不 
同 孔 际 分 布下 , Winkler 和 Pasternak 弹性 支撑 均 可 
提高 模 态 频率 , 且 在 低 环 向 波 数 下 效果 明显 。 这 是 
因为 弹性 支撑 的 作用 相当 于 对 圆柱 壳 施 加 了 约束 ， 
提高 了 结构 刚度 。 相 较 于 Winkler 弹性 支撑 ， 
Pasternak 弹性 支撑 可 较 大 幅度 提高 模 态 频率 。 因 
此 ,在 不 改变 结构 尺寸 时 ,选择 Pasternak 弹性 支撑 
能 调节 微 圆柱 壳 模 态 频率 ,以 避免 结构 发 生 共振 。 


表 3 环 向 波 数 对 无 量 纲 模 态 频率 的 影响 


Tab.3 Effects of circumferential wave number on dimensionless modal frequency 


FGPMs- I 


FGPMs- II 


(K,,K,)=(0,0) (K,,K,)=(1000,0) (K,,K,) 2(1000,1) 


wp wp wp 


(K,,K,) (0,0) (K,,K,) 2(1000,0) (K,,K,) 2(1000,1) 


w »^*p wp wp 


1 0.868 8 0.909 1 2.0203 
2 0.7138 0.7611 2.2404 
3 0.6360 0.6890 2.5809 
4 0.7202 0.768 2 3.0068 
5 0.9611 0.998 1 3.4899 
6 1.3152 1.3428 4.0131 
7 1.7570 1.7778 4.5662 
8 2.276 1 2.2923 5.1439 
9 2.8683 2.8812 5.7432 
10 3.5318 3.5423 6.3627 
12 5.0684 5.0757 7.6604 
15 7.8892 7.8939 9.7568 
18 11.3171 11.3204 12.046 5 
20 13.679 7 13.6798 13.6892 


4.2 外 部 载荷 的 影响 

表 4~ 表 6 分别 描述 了 热 . 力 、 电 3 种 载荷 单独 
变化 下 ,结构 模 态 频率 的 变化 情况 。 表 4 研究 了 温 
度 对 圆柱 壳 无 量 纲 模 态 频率 的 影响 。 结 果 表 明 : 在 


0.8414 0.8694 2.0100 
0.6824 0.715 7 2.2332 
0.604 1 0.6418 2.576 4 
0.696 4 0.729 8 3.0041 
0.947 7 0.9728 3.488 5 
1.3098 1.3283 4.0125 
1.7574 1.7713 4.5665 
2.2812 2.2919 5.1449 
2.8774 2.8859 5.7450 
3.5445 3.5515 6.3653 
5.0884 5.0932 7.6650 
7.9211 7.9243 9.7765 5 
11.3639 11.3660 12.060 6 
13.7020 13.7020 13. 707 7 


圆柱 壳 的 刚度 ,所 以 增 大 温度 梯度 ,圆柱 壳 的 模 态 频 
率 减 小 。 在 Pasternak 弹性 支撑 下 ,温度 对 模 态 频率 
基本 无 影响 。 这 是 由 于 Pasternak 弹性 支撑 较 大 程 
度 的 提高 了 结构 刚度 , 增 大 温度 梯度 不 易 改 变 结 构 
刚度 。 相 比 无 弹性 支撑 ,Winkler 弹性 支撑 下 温度 对 


无 弹性 支撑 和 Winkler 弹性 支撑 下 ,因为 温度 降低 了 


模 态 频率 影响 较 小 。 尤 其 是 Pasternak 弹性 支撑 下 ， 
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模 态 频率 随 温 度 基 本 无 变化 。 因 此 弹性 支撑 会 降低 
温度 对 模 态 频率 的 影响 。 在 无 弹性 支撑 和 Winkler 
弹性 支撑 下 ,孔隙 体积 分 数 越 大 , 微 圆柱 帝 对 温度 越 
敏感 。 说 明 低 孔隙 体积 分 数 可 相对 较 好 的 抵抗 高 温 
度 梯 度 。 比 较 发 现 ,FGPMs- 工 型 孔隙 分 布 对 温度 梯 
度 更 敏感 。 例 如 , 当 e。=0.2 且 温 度 梯度 从 0 增 大 

表 4 
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到 800 开 时 ,孔隙 分 布 为 FGPMs- 工 型 的 模 态 频率 分 
别 减 小 了 0. 916% 和 0. 85496 ; 孔隙 分 布 为 FGPMs- 
I 型 的 微 圆柱 壳 模 态 频 率 分 别 减 小 了 0. 815% 和 
0.769% 。 究 其 原因 ,两 种 孔 阶 下 圆柱 过 对 温度 的 敏 
感性 不 同 主要 体现 在 结构 的 有 效 热 释 电 常数 上 。 
此 ,高 温 环境 下 工作 时 ,可 考虑 FGCPMs-[[ 型 孔隙 分 布 。 


温度 梯度 对 无 量 纲 模 态 频率 的 影响 


Tab.4 Effects of temperature gradient on dimensionless modal frequency 


FGPMs- I 
€o (K, Kp ) 


FGPMs- II 


T,-0K 7,-200K T,-400K T, -600K 7, -800K T, ,=0K 7,-200K T, -400K T, -600K 7, -800K 


io 


0.2 (0,0) 0.8404 0.8385 0.836 6 0.8346 
(1000,0) 0. 866 6 0.864 8 0.8629 0.8611 
(1000,1) 2.008 2 2.0082 2.0082 2.008 2 

0.4 (0,0) 0.8588 0.8563 0.8537 0.8512 
(1000,0) 0.8929 0.8905 0.888 1 0.8856 
(1000,1) 2.0149 2.0149 2.0149 2.0149 

0.6 (0,0) 0. 894 8 0.891 1 0.8874 0.883 7 
(1000,0) 0.9435 0.9400 0.936 5 0.9330 
(1000,1) 2.0283 2.0283 2.0283 2.0283 


表 5 描述 了 轴 向 力 对 圆柱 壳 无 量 纲 模 态 频率 的 


0.8327 0.8343 0.8326 0.8309 0.8292 0.8275 
0.8592 0.8579 0.8562 0.8546 0.8529 0.8513 
2.0082 2.0060 2.0060 2.0060 2.0060 2.0060 
0.848 7 0.8407 0.8388 0.8369 0.8350 0.8331 
0.8832 0. 8669 0.8651 0.8632 0.8614 0.8595 
2.0149 2.008 6 2.008 6 2.008 6 2.008 6 2.008 6 
0.8799 0.848 7 0.846 6 0.8444 0.8423 0.8401 
0.9294 0.8784 0.8763 0.8742 0.8721 0.8701 
2.0283 2.0116 2.0116 2.0116 2.0116 2.0116 


的 刚度 , 轴 癌 压力 减 小 了 圆柱 壳 的 刚度 。Pasternak 


影响 情况 。 当 圆柱 壳 轴 向 受 拉 时 , 模 态 频率 随 轴 向 
力 的 增 大 而 增 大 ; 当 轴 向 受 压 时 , 模 态 频率 随 轴 向 力 
的 增 大 而 减 小 。 这 是 因为 轴 向 拉力 增 大 了 微 圆 柱 壳 


弹性 支撑 下 圆柱 充 对 轴 向 力 敏感 性 最 弱 , Winkler 弹 
性 支撑 次 之 ,无 弹性 支撑 时 人 敏感 性 最 强 。 考 虑 弹性 
支撑 会 降低 轴 向 力 对 模 态 频率 的 影响 。 


表 5 轴 向 力 对 无 量 纲 模 态 频率 的 影响 


Tab.5 Effects of axial load on dimensionless modal frequency 


FGPMs- I FGPMs- Il 
eo (K,,K,) 
N=-1 N=-0.5 N=0 Nm=05 N=1 N=-1 N=-0.5 N=0 N=0.5 N=1 
0.2 (0,0) 0.6630 — 0.7258 0.7837 0.8375 0.8881 0.6775 0.7326 0.7838 0.8318 0.8772 
(1000,0) 0.6959 0.7561 0.8118 0.8639 0.9130 0.7063 0.7593 0.8088 0.8554 0.8996 
(1000,1) 1.9483 1.9686 1.9885 2.0082 2.0277 1.95330 1.9709 1.9886 2.0061 2.0235 
0.4 — (0,0) 0.6179 — 0.7059 0.7840 . 0.8550 0.9206 0.6634 0.7262 0.7840 0.8378 0.8884 
(1000,0) 0.6645 0.7470 0.8212 0.8893 0.9525 0.6963 0.7564 0.8121 0.8642 0.9133 
(1000,1) 1.9349 1.9621 1.9887 2.0150 2.0407 1.9487 1.9689 1.9889 2.0086 2.0280 
0.6 — (0,0) 0.5169 — 0.6646 — 0.7850 0.8893 0.9825 0.6446 0.7179 0.7843 0.8455 0.9026 
(1000,0) 0.5972 0.7289 0.8401 0.9383 1.0271 0.6832 0.7527 0.8163 0.8753 0.9305 
(1000,1) 1.9077 . 1.9489 1.9891 2.0283 2.0665 1.9431 1.9663 1.9891 2.0116 2.0338 


&x 。 FGPMs- 了 [型 孔 际 分 布 对 轴 疝 力 更 敏感 。 例 如 , 当 孔 


据 表 5 可 知 ,孔隙 体积 分 数 越 大 , 轴 向 力 对 模 态 


频率 影响 越 大 。 且 无 论 在 轴 疝 拉力 或 轴 向 压力 下 ， 
大 轴 向 载荷 时 孔 际 体积 分 数 对 模 态 频率 的 影响 较 
大 。 因 此 选择 高 孔隙 体积 分 数 及 大 轴 向 载荷 课 调 节 
模 态 频率 。 比 较 发 现 , FGPMs- 工 型 孔隙 分 布 比 


PETRA e, =0. 2 时 , 当 无 量 纲 轴 向 力 从 -1 增加 
到 1,3 种 弹性 支撑 下 ,孔隙 分布 为 FGPMs- T 79585 pi] 
柱 壳 模 态 频率 分 别 增加 了 33.952% 、31. 197% 和 
4. 075% ;孔隙 分 布 为 FGCPMs- 卫 型 的 圆柱 壳 模 态 频 
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率 分 别 增加 了 29. 476% 、27. 368% 和 3. 610% 。 这 
主要 是 因为 两 种 孔 际 下 有 效 弹 性 常数 和 有 效 质量 密 
度 的 贡献 比 不同 而 导致 敏感 性 不 同 。 因 此 ,车 通过 
轴 向 力 来 调节 微 圆 柱 壳 的 模 态 频率 ,可 选 FGPMs- I 
型 孔隙 分布 。 

表 6 分 析 了 电压 对 圆柱 壳 无 量 纲 模 态 频率 的 影 
啊 。 在 无 弹性 支撑 和 Winkler 弹性 支撑 下 ,施加 正 电 
压 时 模 态 频率 减 小 ,施加 负电 压 时 模 态 频率 增 大 。 
这 是 因为 正 电 压 降 低 了 微 圆 柱 壳 的 刚度 ,而 负电 压 
增 大 了 微 圆柱 壳 的 刚度 。 因 此 ,在 这 两 种 弹性 文 撑 
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压 对 模 态 频率 的 影响 。 在 无 弹性 支撑 和 Winkler 弹 
性 支撑 下 , 减 小 孔隙 体积 分 数 会 增 大 电压 对 模 态 频 
率 的 影响 。 且 在 高 外 加 正 电压 下 ,孔隙 体积 分 数 对 
模 态 频率 的 影响 尤其 显著 。 因 此 ,可 通过 选择 低 孔 
隙 体积 分 数 及 高 外 加 正 电压 参数 来 调节 模 态 频 率 。 
在 无 弹性 支撑 和 Winkler 弹性 支撑 下 ,改变 电压 时 
FGPMs- 丁 型 孔隙 下 模 态 频率 的 变化 较 大 。 例 如 , 当 
eo 2 0. 2 时 , 当 电 压 从 -1 增加 到 1, 孔隙 分 布 为 
FGPMs- [型 的 微 圆柱 这 模 态 频 率 分 别 减少 了 
23. 059% 和 22. 073% ;孔隙 分 布 为 FCPMs- 卫 型 的 微 


下 ,可 通过 改变 外 加 电压 的 方向 和 幅 值 来 调节 模 态 
频率 。 而 在 Pasternak 弹性 支撑 下 ,电压 对 模 态 频率 
基本 无 影响 。 此 时 Pasternak 弹性 支撑 起 主导 作用 ， 
通过 改变 外 加 电压 来 调节 模 态 频率 的 效果 不 明显 。 
与 热 和 轴 向 力 情况 下 相同 ,考虑 弹性 支撑 可 降低 电 


圆柱 壳 模 态 频率 分 别 减 少 了 23.259% 和 22. 370% 。 
这 是 由 于 两 种 孔隙 对 有 效 弹 性 常数 的 影响 不 同 而 导 
致 的 。 因 此 ,可 选择 FCPMs- 卫 型 孔隙 分 布 来 调节 模 


表 6 电压 对 无 量 纲 模 态 频率 的 影响 


Tab.6 Effects of applied voltage on dimensionless modal frequency 


FGPMs- I FGPMs- Il 
ej (K,,K,) 
p=-1 N=-0.5 M=0 M=0.5 N=1 N=-1 N=-0.5 N=0 N=0.5 N=1 
0.2 (0,0) 1.085 — 1.08318 — 0.9716 — 0.9071 0.8375 1.0839 1.0070 0.9665 0.9018 0.8318 
(1000,0) 1.1086 — 1.0531 0.9942 0.9314 0.8639 1.1019 1.0460 0.9868 0.9235 0.8554 
(1000,1) 2.0082 2.0082 2.0082 2.0082 2.0082 2.0061 2.0061 2.0061 2.0061 2.0061 
0.4 (0,0) 1.1026 — 1.0465 — 0.9870 | 0.9035 0.8550 1.0886 — 1.020 0.9718 0.9074 0.8378 
(1000,0) 1.1291 1.0744 1.0167 0.9552 0.8893 1.1088 1.0533 0.9944 0.9317 0.8642 
(1000,1) 2.0150 2.0150 2.0150 2.0150 2.0150 2.0086 2.0086 2.0086 2.0086 2.0086 
0.6 (0,0) 1.1312 1.0761 1.0178 0.9558 0.8893 1.0947 1.0383 0.9785 0.9145 0.8455 
(1000,0) 1.1696 1.1166 1.0607 1.004 0.9383 1.1176 1.0625 1.0042 0.9420 0.8753 
(1000,1) 2.0283 2.0083 2.0283 2.0283 2.0083 2.0116 2.0116 2.0116 2.0116 2.0116 


以 上 分 析 表 明 ,两 种 孔 际 分 布下 微 圆柱 壳 的 模 
态 频率 随 温 度 、 轴 向 力 和 电压 的 变化 趋势 相同 ,只 是 数 
值 上 的 差异 ,可 根据 实际 情况 选择 何 种 孔隙 分 布 。 为 
简化 分 析 , 在 以 下 讨论 中 只 考虑 FGPMs-] 型 孔 际 分 布 。 


4.3 不 同 弹性 支撑 和 材料 体积 分 数 的 影响 


图 3 ~5 分 别 为 在 不 同 弹性 支撑 和 材料 体积 4 
数 指数 下 无 量 纲 模 态 频 率 随 温 度 、 轴 向 力 和 电压 的 
变化 曲线 。 

结果 表明 ,无 弹性 支撑 下 ,材料 体积 分 数 指数 越 
大 , 微 圆 柱 壳 的 模 态 频率 越 大 且 随 温度 、 轴 向 力 和 电 
压 的 变化 率 越 小 。 这 是 因为 材料 体积 分 数 指数 越 
大 ,PZT-4 的 占 比 越 大 ,圆柱 完 的 刚度 增 大 。 说 明 增 
X PZT-4 的 占 比 ,圆柱 壳 可 承受 高 温度 ,大 轴 向 载荷 


态 频 率 随 轴 向 力 变化 明显 。 而 改变 温度 梯度 或 外 加 
电压 , 模 态 频 率 基 本 不 变 。 此 时 可 通过 改变 材料 体 
积分 数 指数 来 调节 微 圆柱 壳 的 模 态 频率 ,并 可 以 达 
到 良好 的 效果 。 


4.4 孔隙 体积 分 数 的 影响 


图 6 为 无 量 纲 模 态 频率 随 孔 际 体积 分 数 的 变化 
曲线 。 结 果 表 明 , 当 孔 际 体积 分 数 小 于 0.2 时 , 模 态 
频率 变化 不 大 ; 当 孔 际 体 积分 数 大 于 0.5 时 , 增 大 孔 
际 体 积分 数 模 态 频率 的 变化 趋势 明显 不 同 。 无 论 在 
无 弹性 支撑 下 还 是 在 Pasternak 弹性 支撑 下 , 当 材 料 
体积 分 数 指数 较 小 时 , 模 态 频率 随 孔 际 体积 分 数 的 
增 大 而 逐渐 减 小 ; 当 材 料 体积 分 数 指数 较 大 时 , 模 态 
频率 随 孔 际 体 积分 数 的 增 大 而 增 大 。 究 其 原因 ,这 


以 及 高 电压 。 在 Pasternak 弹性 支撑 下 ,圆柱 壳 的 模 


可 以 从 孔隙 体积 分 数 改 变 有 效 弹性 模 量 和 有 效 密度 
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图 3 无 量 纲 模 态 频率 随 温度 梯度 的 变化 曲线 


Fig.3 Vibration of dimensionless modal frequency 


with temperature gradient 
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5 无 量 纲 模 态 频率 随 电 压 的 变化 曲线 


Fig.5 Vibration of dimensionless modal frequency 
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with applied voltage 
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的 角度 来 说 明 。 当 材料 体积 分 数 指数 较 小 时 , 微 圆 
柱 壳 的 刚度 较 小 。 此 时 增 大 孔 际 体积 分 数 ,造成 对 
有 效 弹性 模 量 的 削减 大 于 有 效 密度 ; 当 材 料 体积 4 
数 指数 较 大 时 , 增 大 孔 际 体积 分 数 ,显著 降低 了 微 圆 


3.0 
一 一 (Ko K,)-(1 000, 1), N-0.1 


35b —9-(K.. K)-( 000, 1), N-1 


2.0 


—o— (K,, K;)-(1 000, 1), N-inf 
—a— (K,, K,)-(0, 0), N=0.1 


无 量 纲 模 态 频率 Q 


05| a (K, K)-(0,0),N-1 
—e— (K,, K,)=(0, 0), N-inf 


0 01 02 03 04 05 06 07 03 
孔隙 体积 分 数 e 


图 6 无 量 纲 模 态 频率 随 孔隙 体积 分 数 的 变化 曲线 
Fig.6 Vibration of dimensionless modal frequency 


with porosity volume fraction 


4.5 ”过 体 厚度 的 影响 


图 7 和 图 8 分 别 展示 了 基于 经 典 连续 体 理论 和 


修正 的 偶 应 力 理论 下 , 微 圆柱 壳 无 量 纲 模 态 频率 随 
厚度 参数 的 变化 情况 。 
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图 7 无 量 纲 模 态 频率 随 壳 体 厚度 的 变化 曲线 (1 =0) 


Fig.7 Vibration of dimensionless modal frequency 


with shell thickness (7-0) 
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Fig.8 Vibration of dimensionless modal frequency 


with shell thickness (/ 23A) 
结果 表明 ,基于 修正 的 偶 应 力 理论 的 模 态 频率 
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大 于 基于 经 典 连 续 体 理论 的 模 态 频率 。 说 明 考 虑 材 
料 特 征 长 度 会 使 得 微 圆 柱 壳 的 刚度 变 大 。 随 着 充 体 
厚度 增加 ,有 无 弹性 支撑 之 间 的 模 态 频率 逐渐 减 小 ， 
TE 0.01 « h/R <0. 02 时 , 微 圆柱 壳 的 模 态 频率 差异 
明显 。 说 明 Pasternak 弹性 文 撑 对 薄 壳 的 影响 更 大 。 
究 其 原因 , 薄 壳 较 “ 软 ”, 弹性 文 撑 对 其 作用 效果 显 
著 。 对 于 不 同 长 度 的 微 圆 柱 克 ,对 比 其 在 有 无 弹性 
支撑 下 的 模 态 频率 ,在 1<LR<5 之 间 , 发现 
Pasternak 弹性 支撑 提高 短 圆柱 过 的 模 态 频率 更 加 显 
著 。 这 主要 是 因为 Pasternak 弹性 文 撑 中 的 剪 切 层 
对 短 圆 柱 壳 的 作用 效果 更 明显 。 


5 f idt 


基于 修正 的 偶 应 力 理论 和 三 阶 剪 切 变 形 壳 体 理 
论 推导 了 弹性 支撑 功能 梯度 压 电 多 孔 微 圆柱 壳 的 运 
动 方程 ,分 析 了 各 种 参数 对 微 圆 柱 壳 模 态 频率 的 影 
响 。 主 要 结论 如 下 。 

1) Winkler 弹性 支撑 仅 影 响 微 圆 柱 壳 的 径 向 振 
动 ,而 Pasternak 弹性 支撑 影响 径 回 和 环 向 振动 。 与 
Winkler 弹性 支撑 相 比 ,Pasternak 弹性 支撑 可 较 大 幅 
度 提高 微 圆 柱 壳 的 模 态 频率 。 

2) 改变 弹性 支撑 的 刚度 系数 、 轴 向 力 ` 外 加 电 
压 、 和 孔隙 分 布 类 型 .孔隙 体积 分 数 、 材 料 体积 分 数 指 
数 和 结构 尺寸 等 , 均 可 起 到 调节 微 圆柱 壳 模 态 频率 
的 作用 。 温 度 . 轴 向 压力 和 外 加 正 电 压 均 降低 微 圆 
柱 壳 的 刚度 ,而 轴 向 拉力 和 外 加 负电 压 增 大 微 圆柱 
壳 的 刚度 。 

3) 增 大 和 孔 际 体积 分 数 ,温度 载 答 或 外 加 轴 向 力 
对 模 态 频率 的 影响 增 大 ,外 加 电压 对 模 态 频率 的 影 
响 减 小 。 在 不 同 的 材料 体积 分 数 指数 下 , 增 大 孔 院 
体积 分 数 模 态 频 率 的 变化 趋势 不 同 。 

4) 弹性 支撑 会 降低 温度 、 轴 向 力 和 电压 对 模 态 
频率 的 影响 , 对 于 Pasternak 弹性 支撑 而 言 尤 其 明 
o Pasternak 弹性 支撑 对 短 圆柱 过 或 薄 圆 柱 壳 的 模 
态 频率 影响 更 加 显著 。 
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